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Abstract  — To  truly  optimize  the  deployment  of  DoD 
assets,  there  exists  a  fundamental  need  for  predictive 
tools  that  can  reliably  estimate  the  current  and  reasonably 
predict  the  future  capacity  of  complex  systems. 
Prognosis,  as  in  all  true  predictions,  has  inherent 
uncertainty,  which  has  been  treated  through  probabilistic 
modeling  approaches.  The  novelty  in  the  current 
prognostic  tool  development  is  that  predictions  are  made 
through  the  fusion  of  stochastic  physics  -of-failure 
models,  relevant  system  or  component  level  health 
monitoring  data  and  various  inspection  results. 
Regardless  of  the  fidelity  of  a  prognostic  model  or  the 
quantity  and  quality  of  the  seeded  fault  or  run -to -failure 
data,  these  models  should  be  adaptable  based  on  system 
health  features  such  as  vibration,  temperature,  and  oil 
analysis.  The  inherent  uncertainties  and  variability  in 
material  capacity  and  localized  environmental  conditions, 
as  well  as  the  realization  that  complex  physics -of-failure 
understanding  will  always  possess  some  uncertainty,  all 
contribute  to  the  stochastic  nature  of  prognostic 
modeling.  However,  accuracy  can  be  improved  by 
creating  a  prognostic  architecture  instilled  with  the  ability 
to  account  for  unexpected  damage  events,  fuse  with 
diagnostic  results,  and  statistically  calibrate  predictions 
based  on  inspection  information  and  real-time  system  level 
features.  In  this  paper,  the  aforementioned  process  is 
discussed  and  implemented  first  on  controlled  failures  of 
single  spur  gear  teeth  and  then  on  a  helical  gear  contained 
within  a  drivetrain  system.  The  stochastic,  physics -of- 
failure  models  developed  are  validated  with  transitional 
run-to -failure  data  developed  at  Penn  State  ARL.  Future 
work  involves  applying  the  advanced  prognostics  process 
to  helicopter  gearboxes. 

Table  of  Contents 

1.  Introduction 

2.  RUN-T  O-FAILURE/SEEDED  FAULT 


Testing 

3.  Gear  Prognostics 

4.  Results 

5.  knowledge  Fusion 

6.  Conclusions 

1.  Introduction 

Battlefield  commanders  currently  have  limited  information  to 
evaluate  their  assets’  operational  capability  or  the  next 
mission’s  effect  on  that  capability.  Thus,  they  must  make 
“go”/“no  go”  decisions  that  can  result  in  the  underutilization 
of  assets  with  some  remaining  capability  and  the  risk  of 
overtaxing  deployed  assets  that  have  insufficient  capability  to 
meet  mission  requirements.  The  assessment  of  current  health 
state  and  the  prediction  of  future  capability  for  complex 
mechanical  systems  is  a  core  enabling  technology  for  such 
decision  aiding. 

The  focus  of  this  current  effort  was  to  demonstrate  the 
methodology  and  develop  tools  to  provide  state  awareness 
and  capability  prediction  on  a  mechanical  system.  The 
approach  follows  an  integrated  mathematical  (probabilistic) 
framework  for  utilizing  ALL  of  the  available  data/information 
that  can  lead  to  better  predictions  (less  uncertainty)  on  future 
capability  of  a  component/machine.  This  includes  material 
level  fatigue  models,  system  level  feature  models  and  raw 
health  monitoring  measurements  that  provide  the 
data/interaction  between  system  and  component  levels.  This 
mathematical  framework  requires  a  consistent  normalization 
and  probabilistic  updating  process/algorithm  that  occurs  at 
specified  triggering  points  in  the  life  of  a  component  based  on 
indication  from  features.  These  features,  extracted  from  various 
system  level  measurements,  can  be  used  to  understand  the 
progression  of  a  failure  (potentially  with  competing  failure 
modes).  Considerable  work  and  useful  results  have  been 
identified  in  this  technical  area,  primarily  from  the  CBM  and 
machinery  health  monitoring  communities  and  the  intent  of  this 
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paper  is  not  to  reproduce  those  results.  Rather,  the  focus 
of  this  paper  is  on  maximizing  the  utility  of  features  based 
on  their  ability  to  track  material  damage  progression 
through  particular  failure  modes  and  provide  statistical 
correlation  with  actual  component  damage  levels  that  can 
then  be  linked  with  the  failure  mode  material  damage 
models. 


2.  Run-to- Failure/Seeded  Fault  Testing 


Developing  validated  models  for  predicting  asset 
readiness  requires  normal,  seeded,  and  transitional  (run- 
to-failure)  data  that  have  well-documented  ground  truth 
information.  In  general,  there  is  a  lack  of  existing,  high 
fidelity  data  sets  pertaining  to  fault  initiation  and 
evolution.  Thus  the  desire  for  accurate  prognostic  tools 
drives  the  accompanying  need  for  realistic  data  sets  to 
develop  asset  state  prediction  technology  (e.g.,  failure 
models,  diagnostic  algorithms,  predictive  techniques,  data 
fusion,  reasoning). 

)  As  illustrated  in  Figure  1,  Run-To -Failure/Seeded  Fault 
(RTF/SF)  results  are  used  to  1)  provide  realistic 
indications  of  state  awareness,  that  is,  health  indicators 
(i.e.  vibration  features)  and  non-destructive  evaluations 
related  to  actual  material  damage  levels  2)  identify 
sensitivities  and  uncertainties  in  the  effects  of  material 
properties  and  manufacturing  defects  on  component 
capacity  and  finally,  3)  provide  invaluable  calibration  of 
any  range  of  prognostic  model  at  various  times  in  the  life 
of  a  component  so  that  it  can  be  evaluated  both  in  terms 
of  long-term  capability  prediction  (asset  management)  and 
more  near-term  damage  minimization  (fault 
accommodation). 


Figure  1  -  Capability  Prediction 

3.  Gear  Prognostics 


application,  the  purpose  of  the  effort  was  to  very  accurately 
correlate  a  2-D  Finite  Element  fracture  mechanics  model  and 
associated  crack  initiation  and  propagation  lifing  algorithms  to 
tooth  stiffness  and  acoustic  emissions  changes.  This  is 
achieved  by  updating/adapting  material  property  distributions 
or  choice  of  lifing  algorithm  during  damage  progression  (used 
in  the  local  strain  stochastic  fatigue  models)  based  on 
measured/inferred  conditions  at  the  local  failure 
location/region.  The  second  prognostic  model  was  developed 
for  a  helical  gear  in  an  industrial  gearbox  run  to  failure  at  ARL. 
Using  a  high  fidelity  model  (build  upon  contact  element  and  3- 
D  fracture  mechanics  FE  models  of  the  gear  and  a  6  DOF 
torsional  model  of  the  drivetrain)  and  well-proven  vibration 
features,  an  attempt  was  made  at  predicting  current  and  future 
material  level  damage  (capability)  as  a  function  of  torque  from 
system  level  vibration.  Thus,  system  level  features  and 
operational  state  awareness  are  used  to  ascertain  current  and 
future  material  capacity. 


Single  Spur  Gear  Tooth 


Figure  2  shows  a  schematic  of  the  gear  tooth  testing  machine 
used  to  fail  spur  gear  teeth. 


Figure  2  -  Single  Tooth  Test  Rig 

In  this  study,  the  material  was  ASM -6265  and  the  teeth  were 
case  hardened  to  approximately  60-62  HRC  to  a  depth  of  0.036 
-  0.044  inches.  During  the  fatigue  tests,  tooth  loading  was 
along  a  line  of  contact  at  the  highest  point  of  single  tooth 
contact  and  the  gear  tooth  was  cyclically  loaded  using  a 
Haversign  profile  with  alternating  loads  ranging  from  300-800 
min.  to  7,200-11,000  max.  Each  test  was  faulted  prior  to 
complete  tooth  failure  based  on  a  tooth  deflection  limit.  An 
ANSYS  model  of  the  tooth  was  used  to  calculate  true  strain  at 
the  tooth  root  for  a  number  of  different  test  cases.  As  shown 
in  Figure  6,  each  test  case  was  plotted  on  a  strain -life  curve  to 
empirically  determine  the  fatigue  parameters  n,c,K,and  £f  used 
in  the  Eq.  #1.  The  banding  shown  illustrates  approximated  +/- 
3a bounds  on  the  low,  high,  and  combined  S-N  curves. 

1  \—]  [— 1  f-1  (!) 

NflL  =  — -[G  L  (true)- G  (mean)  p(n  c)^  •  K^-(n  c)J  •  Ef^ c ' 


Under  this  program,  two  gear-based  prognostic  models 
were  developed.  The  first  was  for  a  single  spur  gear  tooth 
and  validated  Run-to-Failure  results  and  deflection  and 
acoustic  emission  data  from  a  simple  test  rig  at  Penn  State 
ARL.  While  quite  removed  from  a  realistic  gearbox 


NflL  =  the  LCF  life  for  the  gear  (L) 

aL(true)  =  localized  true  plastic  stress  amp.  at  a  tooth  root 
n  =  cyclic  strain  hardening  exponent 
c  =  fatigue  ductility  exponent 
K  =  cyclic  strength  coefficient 


Ef=  fatigue  ductility  coefficient 
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Figure  3  -  S-N  curve  with  mat.  property  uncertainties 

This  model  was  run  as  a  Monte  Carlo  simulation  whereby 
random  samples  were  taken  from  each  of  the  material 
property  distributions  and  then  randomly  combined  in  the 
above  equation.  The  results  of  this  statistical  crack 
initiation  model  for  a  mean  loading  of  4950  lbf.  is  shown  in 
Figure  4.  The  actual  crack  initiation  occurred  around 
3.32E8  total  accumulated  lbs.  (as  inferred  from  the  acoustic 
emission  level  and  tooth  deflection  shown  in  Figure  5)  or 
approx.  70K  cycles/140K  reversals. 


Figure  4  -  PDF  on  crack  initiation 
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Figure  5  -  Tooth  Deflection  vs.  Accumulated  Load 

The  FE  model  automatically  propagated  a  crack  at  fixed 
increments  perpendicular  to  the  maximum  principle  stress  if  an 
initial  notch  was  created.  This  model  provided  the  orthogonal 
threshold  intensity  ranges  at  the  crack  front  for  the  Foreman 
model  given  in  Eq.  2. 


Figure  6  -  2-D  Finite  Element  Model 

The  Forman  equation  states  that  for  Phase  II  and  Phase  III 
crack  growth  with  mean  stress  effects: 


da  _  CAKm  (2) 

~dN~(\-R)K:ril-AK 


A  fracture  mechanics  stochastic  model  and  associated  2-D 
FE  model  was  also  created  for  this  gear  (Figure  6).  While 
the  cycles  from  “like  new”  condition  to  crack  initiation  is 
much  less  that  the  cycles  from  crack  initiation  to  unstable 
crack  growth,  the  purpose  of  this  effort  was  to  show  that  a 
life  prediction  can  be  updated  based  on  state  awareness 
(crack  detected  or  core -to -case  material  transition). 


where: 

a  =  crack  length 
N  =  Number  of  cycles 
AK  =  Threshold  intensity  range 
C,m  =  Phase  II  fracture  characteristics 
R  =  mean/alternating  stress  ratio 

The  threshold  intensity  range  for  a  given  crack  length,  coupled 
with  gross  estimates  of  mean  “C”  and  “m”  value  of  6.86E-9  and 
2.8  respectively,  was  used  to  obtain  distributions  on  the 
number  of  stress  cycles  to  the  crack  length  at  which  the  test 
was  halted  with  a  crack  of  0.083”  in  81,820  cycles.  Thus  it  can 
be  estimated  that  the  crack  propagated  from  microscopic  size 
to  0.083  inches  in  approx.  24K  reversals,  which  also  falls  within 


the  model-based  results  for  crack  propagation  (Figure  7). 


Figure  7  -  Crack  propagtion  PDF 

In  this  model,  the  total  probability  of  failure  is  defined  as 
the  Probability  of  crack  initiation  (Pi)  *  the  Probability  of 
the  crack  propagating  to  unstable  crack  growth  (Pp).  The 
observed  features  clearlymarked  the  transition  from  crack 
initiation  to  crack  growth  and  provided  a  critical  piece  of 
information  for  the  prognostic  model.  Specifically,  when 
the  system  features  suggest  a  crack  has  initiated  the 
probability  of  crack  initiation  can  be  estimated  as  close  to 
unity,  therefore  reducing  the  uncertainty  in  the  total 
probability  of  failure  prediction  significantly.  In  addition, 
the  predicted  stress  intensity  factor  used  in  the  model  can 
be  updated  based  on  the  relationship  between  the 
measured  deflection  and  inferred  crack  length/location 
prediction.  This  update,  based  on  the  inferred  knowledge 
on  crack  length/location,  will  allow  the  model  to  remain 
“informed”  on  the  current  state  of  damage  so  that 
improved  predictions  on  remaining  capability  can  be  made. 

In  the  case  study  developed  thus  far,  immediately  prior  to 
the  detected  crack  initiation,  the  total  probability  of  failure 
could  only  be  calculated  as  Pi  (15%)  *  Pp(99.43%)  = 
14.91%  .  As  soon  as  the  crack  is  detected,  the  Probability 
of  crack  initiation  (Pi)  is  assigned  a  value  of  unity  and  the 
Probability  of  failure  due  to  propagation  (Pp)  is  reset  to 
zero.  Thus  the  uncertainty  in  the  current  component 
capability  is  reduced  by  nearly  an  order  of  magnitude. 

Finally,  a  demonstration  was  developed  to  illustrate  how 
crack  initiation  and  propagation  and  case  to  core 
propagation  transition  can  influence  the  current  and  future 
failure  probabilities  (Figure  8). 

In  this  demo,  three  outputs  are  given  to  evaluate  the  crack 
progression:  the  current  total  probability  of  failure 
obtained  from  the  current  cumulative  distribution  function, 
the  future  probability  based  on  the  entered  projected 
cycles  and  the  current  remaining  number  of  cycles  to 
failure  with  a  95%  confidence  interval. 


Figure  8  -Single  Tooth  Prognostics  demo 
Gearbox  application  -  Helical  gear 

Another  prognostic  model  was  developed  for  a  helical  gear  in 
an  industrial  gearbox.  Details  on  the  test  rig  used  to  generate 
Run-to -Failure  data  may  be  found  in  [3].  The  prognostic 
concept  employed  here  is  similar  to  the  single  gear  tooth  model 
already  described  in  the  sense  that  system  level  features  are 
used  to  improve  the  accuracy  of  component  level  lifing  models, 
but  different  in  terms  of  level  of  modeling  fidelity  and  relevance 
to  real-world  failures. 

Figure  9  outlines  the  architecture  of  this  prognostic  model. 
Standalone  physics-of-failure  models  predict  current  and 
future  material  damage  levels  as  a  function  of  speed  and 
torque.  If  available,  diagnostic  results  from  system  level 
features  (i.e.  crack  detected)  and  NDE  inspection  results  are 
used  in  a  probabilistic  fusion  process  to  update/calibrate  the 
model  predictions. 


Figure  9  -  Flowchart  for  Helical  Gear  model 

The  encapsulated  prognostic  models  were  built  based  on  3 
different  FE  models  and  2  stochastic  lifing  algorithms.  The 
models  were  as  follows: 

1.  2-D  FE  model  w/contact  elements 

2.  Purpose:  Determine  the  torsional  forcing  function  as  a 


function  of  rotation  and  effective  crack  length. 

3.  6  DOF  Torsional  model 

4.  Purpose:  Transform  the  torsional  forcing  function 
into  simulated  axial  vibration. 

5.  3-D  FE  of  a  cracked  gear  tooth 

6.  Purpose:  Obtain  equivalent  2-D  gear  stiffness 

reduction  as  a  function  crack  length  and  stress 
intensity  ranges  at  the  crack  front  (Figure  10). 


Figure  10  -3-D  modeling  of  a  cracked  tooth 


While  quite  advanced,  several  assumptions  were  still 
made  in  the  modeling  effort.  Principally,  the  3-D  crack  was 
assumed  planar  which  was  deemed  a  reasonable 
assumption  for  phase  II  crack  growth  the  experience  with 
the  failures.  Secondly,  the  contact  element  analysis  was  2- 
D  and  thus  did  not  accurately  reflect  the  forcing  function 
created  from  a  helical  gear  set  in  mesh  but  still  provided  an 
improvement  in  the  calculated  forcing  function  over  basic 
sinusoidal  forcing  functions. 

Figure  1 1  attempts  to  illustrate  the  interaction  between  the 
models  both  for  the  development  of  the  encapsulated 
prognostic  model  and  when  it  is  functioning.  It  is 
interesting  to  note  that  vhen  the  model  is  in  use,  it 
produces  estimates  of  material  damage  (life  used  /  life 
remaining)  for  a  given  mission  profile  (speed,  torque). 
These  predictions  are  updated  based  on  how  the 
observed  vibration  response  differs  from  the  expected 
response  and  how  that  difference  relates  to  run-to -failure 
experiences. 
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Figure  11  -  Model  process  development  and  implementation 

flow 

4.  Results 

Fundamentally,  the  prognostic  model  produces  a  total 
probability  density  function  (PDF)  as  a  function  of  time/cycles 
based  on  probability  of  initiation  times  the  probability  of 
propagation.  This  distribution  with  two  actual  failures 
superimposed  is  given  in  Figure  12.  In  application,  the  model 
would  yield  a  97%  confidence  interval  on  the  time/cycles 
remaining.  One  test  shown  failed  nearly  at  this  interval. 


Figure  12  -  Validation  on  2  runs  at  300%  load 


Figure  13  shows  a  plot  of  the  damage  as  a  function  of  time. 
The  bars  illustrate  97%  confidence  bounds  on  the  damage 
level  at  specified  times.  One  recorded  failure  occurred  37  hours 
after  a  300%  load  was  applied.  The  total  elapsed  time  including 
a  90-hr.  break  in  period  was  127  hrs.  From  Figure  13,  the 
probability  of  tooth  failure  was  51%  at  127  hrs. 


Figure  13  -  Damage  level  prediction  vs.  time 

5.  Knowledge  Fusion 
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And:  m(A)  +  m(B)  {True}  =  (0.13  +  0.13  +0.06)/(l-(0.01 

+0.39))  =  0.53 

This  result  is  called  the  “belief’  and  it  is  the  fused  probability 
lower  bound.  The  uncertainty  in  this  result  is  the  following: 

m(A,A’)  +  m(B,B’)  =  0.06  /  (1-(0.01  +  0.39))  =  0.10  or  +/-  0.05 

Hence,  the  probability  of  Fault  A  having  actually  occurred 
given  the  diagnostic  output  and  in-field  experience  is  0.58  +/- 
0.05. 


While  the  details  on  the  probabilistic  fusion  process 
alluded  to  in  Figure  9,  are  outside  the  scope  of  this  paper 
(more  detail  can  be  found  in  [7]),  fusion  with  the 
Dempster-Shafer  method  is  provided  here. 


In  the  Dempster-Shafer  approach,  uncertainty  in  the 
conditional  probabilities  is  considered.  The  Dempster- 
Shafer  methodology  hinges  on  the  construction  of  a  set, 
called  the  frame  of  discernment,  which  contains  every 
possible  hypothesis.  Every  hypothesis  has  a  belief 
denoted  by  a  mass  probability  (m).  Beliefs  are  combined 
in  the  following  manner. 


Belief  (Hn) 


AnB=Hn _ 

1-  -rtij  (B) 

AnB= 0 


(3) 


The  technique  can  be  best  explained  through  the  use  of 
the  following  example. 


Dempster-Shafer  fusion  as  well  as  Baysian  combination  and  a 
simple  weigthed  a/erage  were  performed  on  seven  vibration 
features  out  of  26  from  a  run-to-failure  test.  As  shown  Figure 
14,  the  features  indicated  a  crack  at  108  hrs  (index  269  on  the 
graph)  with  high  confidence  over  a  very  short  time. 


Figure  14  -  Feature  Fusion 


Given: 

A  diagnostic  classifier  detects  Fault  A  with  the  following 
probability  and  associated  uncertainty: 

PA  =  0.80+/-  0.15 

The  a  priori  probability  of  Fault  A  occurring  (based  on 
current  conditions  and  a  priori  information)  is  the 
following: 

PB  =  0.30  +/-  0.10 

Therefore:  m(A)  =  0.65  m(A’)  =  0.05 

m(A,A’)  =  0.30 

m(B)  =  0.20  m(B’)  =  0.60 

m(B,B’)  =  0.20 


The  current  probability  of  crack  initiation  is  a  calculated  from 
the  fusion  of  features  with  the  model.  Thus  when  the  fused 
features  alarm  that  a  crack  has  occurred  with  high  confidence, 
the  current  probability  of  failure  is  just  the  probability  of 
propagation  to  failure.  This  knowledge  fusion  allows  the 
uncertainty  in  the  prognostic  model  to  be  greatly  reduced. 

6.  Conclusions 

The  authors  demonstrated  a  methodology  and  prognostic  tool 
development  approach  t>  provide  accurate  state  awareness 
and  improved  capability  prediction  for  machinery  systems. 
Two  examples  of  prognostic  models  have  been  shown  based 
on  a  prognostic  architecture  that  allows  features  to 
update/calibrate  stochastic  physics-of-failure  models.  Using 
system  level  observable  features  to  orient  and  update  the 
models  tracking  material  condition  provides  a  sound  approach 


for  improved  awareness  of  current  health  state.  The  real¬ 
time  calibrated  model  is  then  available  to  predict  future 
performance  over  a  range  of  potential  loads  and 
environmental  conditions.  The  development  of 
prognostic  tools  such  as  those  described  here  can  be 
accomplished  through  the  application  of  material  damage 
models  to  run -to -failure  data  sets.  Deployment  of 
calibrated  prognostic  tools  within  a  tactical  decision- 
aiding  environment  will  provide  intelligent  asset 
management  and  lead  to  improved  readiness  and  mission 
assurance. 
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